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When tensioned, ordinary materials expand along the direction of the applied force. Here, 
we explore network concepts to design metamaterials exhibiting negative compressibility 
transitions, during which a material undergoes contraction when tensioned (or expansion 
when pressured). Continuous contraction of a material in the same direction of an applied 
tension, and in response to this tension, is inherently unstable. The conceptually similar 
effect we demonstrate can be achieved, however, through destabilisations of (meta)stable 
equilibria of the constituents. These destabilisations give rise to a stress-induced solid-solid 
phase transition associated with a twisted hysteresis curve for the stress-strain relationship. 
The strain-driven counterpart of negative compressibility transitions is a force amplification 
phenomenon, where an increase in deformation induces a discontinuous increase in response 
force. We suggest that the proposed materials could be useful for the design of actuators, 
force amplifiers, micro-mechanical controls, and protective devices. 

Metamaterials, engineered materials that gain their properties from structure rather than com- 
position, have led to the study of a myriad of properties not exhibited by most (and in some cases 
all) natural materials. These new properties are typically characterised by negative constants, mod- 
uli, or indexes. In the case of electromagnetic properties, a revolution started a decade ago with 
the experimental realisation of Veselago's predictions [lj on hypothetical materials with negative 
refractive index El[3l, and related applications such as to the design of sub- wavelength lens [4] 
and electromagnetic cloaks [ 5 J. Similar concepts have also been explored to create materials with 
unusual acoustic properties, including negative effective mass density for acoustic focusing (6). 
These properties contrast with those of most conventional materials. 

Conventional materials are also known to deform along the direction of an applied force in such 
a way that they expand when the force is tensional and contract when it is compressive, regard- 
less of how they respond in transverse directions. Thus, one can ask whether a material could be 
designed to exhibit the opposite behaviour, corresponding to negative compressibility. Previous 
studies have shown that porous materials can expand under hydrostatic pressure by sorption or 
infiltration of particles of the medium 117141111 and that resonant composites can exhibit negative 
effective elastic constants in response to acoustic waves [IT2l - [T4ll . Similarly, it has been shown 
that one component of a composite material (but not the material itself) can exhibit decreasing re- 
sponse force as a function of increasing applied deformation lfT5ll . These negative responses have 
led to important applications, such as composites with unusually large stiffness lfT6l or acoustic 
shielding potential [12], but they rely on the system being either mechanically open or dynami- 
cally excited. The question being raised here, however, is whether the material itself (a thermody- 
namically closed system) can exhibit decreasing response deformation as a function of increasing 
applied force (changed over macroscopic timescales). 

The fundamental problem with this question is that closed systems are thermodynamically 
forbidden from expressing such negative compressibility ifTTl . Indeed, for a closed system in 
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equilibrium, if increase in force leads to deformation opposing the applied force, then the origi- 
nal configuration is necessarily unstable. This has been interpreted as an indication that negative 
compressibility in the direction of the applied forces cannot be realised for any closed natural or 
engineered material. However, this argument tacitly assumes that the changes in force are small 
and that the equilibrium survives this change. While true for infinitesimal changes, here we show 
that this is not necessarily true if the force is changed by a finite amount, as expected in many real- 
istic situations. In this case, the equilibrium can become unstable or disappear, forcing the system 
to move to a different equilibrium. We exploit this mechanism to design and theoretically analyze 
materials that have this property, and for which the new equilibrium corresponds to a contraction 
(dilation) when the system is tensioned (pressured), as illustrated in Fig. [T^,b. We refer to this 
discontinuous change as a negative compressibility transition. This transition manifests itself as 
an infinite negative value of the differential compressibility, which is also suitably expressed as a 
finite negative value of the finite-difference compressibility (see Table SI for a summary of def- 
initions). Note that this is a bulk property, since volume can decrease (increase) under increased 
tension (pressure) due to the material's negative response along each and every direction of the 
applied external force. 

Negative compressibility transitions are distinct from negative Poisson's ratio HHEl and neg- 
ative normal stress |[20l . where the unusual response is transverse to the applied force. They are 
also distinct from negative incremental stiffness HUE), characterised by a decrease in the result- 
ing restoration force for increasing deformation. Finally, the proposed negative compressibility 
should be contrasted with stretch densification and the related phenomena of linear and area ex- 
pansion induced by hydrostatic pressure in closed systems J71[22l|23]], which are compensated by 
contraction in other directions, leading to ordinary behaviour for the bulk modulus. Figure [ij-e 
contrasts negative compressibility transitions with negative Poisson's ratio, negative incremental 
stiffness, and stretch densification. 

In order to demonstrate that materials can exhibit negative compressibility transitions, we con- 
sider constituent elements formed by a system of four particles separated by distances x, y, z and 
h and interacting via general potentials V X9 V y , V z and 14, as indicated in Fig. [2^. This system is 
characterised by the total potential 

V(x, y, h, F) = V x (x) + V y (y) + V z (y-x) + V x (h - y) + V y (h - x) + V h (h) - Fh, (1) 

where F is an externally applied force (positive for tension and negative for pressure). We also 
assume there is energy dissipation determined by contact with a heat bath, so that the system can 
reach equilibrium after being perturbed. As in the case of all previously studied closed systems, 
negative compressibility cannot occur continuously (Methods Summary). However, negative com- 
pressibility transitions are still possible in this system by exploiting a bifurcation similar to the one 
that occurs for the potential U{^T) = — £ 3 /3+ J 7 ^, where £ is a variable and T an external param- 
eter. For T > 0, this potential has a stable equilibrium point at = — \[T and an unstable one at 

= \[T\ for T = there is a single degenerate equilibrium point at = 0; for T < there is no 
equilibrium point. Therefore, if T is decreased from positive through zero, the stable equilibrium 
point vanishes giving rise to a destabilisation, and the system responds discontinuously. 

In the case of our system, it is relevant to destabilise the V z bond by increasing (decreasing) the 
force F so as to drive the system to a different equilibrium corresponding to smaller (larger) h. 
This can be shown to be possible for potentials satisfying suitable conditions. We assume F to be 
a tensional force, as reciprocal results can be achieved for pressures. 

The destabilisation weakens the V z bond and causes the system to transition from a predomi- 
nantly series V x — V z — V x configuration (which we call the coupled configuration) to a predomi- 
nantly parallel V x — V y , V y — V x configuration (the decoupled configuration). This causes the bonds 
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defined by V x to contract as the stress becomes distributed between them. If the shortening of x 
overcomes the lengthening of y, the system as a whole contracts (i.e., h decreases) in response to 
increasing tension (Fig. [2^). This is achieved if the potential V x is "soft" in the neighbourhood of 
the transition but "hardened" for smaller x and/or the potential V y is "soft" in the neighbourhood 
of transition but "hardened" for larger y. This both illustrates that harmonic potentials cannot lead 
to the desired behaviour in such one-dimensional systems (although they can in two-dimensional 
configurations — see Supplementary Information, Fig. S2) and provides a criterion for the selection 
of the potentials (see Methods Summary). 

Our design of this system is inspired by the key observation that the realised equilibrium is 
not necessarily optimal in a decentralised network. The best known precedent to this is the in- 
sight from Braess G4ll that adding a road to a traffic network may increase rather than decrease 
the average travel time. Related concepts have been considered in connection with computer sci- 
ence and game theoretical studies E5L and have long been part of the transportation [26] and 
economics ll27l literature, amounting to over 90 years of multidisciplinary research. Analogous 
effects have also been identified in physical networks, including increase of current upon the re- 
moval of an intermediate conductor in electric networks and increase of restoring force upon the 
removal of a support string in spring-string networks [28] . All these are examples in which the 
equilibrium realised by the system can be brought closer to the optimum by constraining the struc- 
ture of the network. Our system is devised such that a conceptually similar phenomenon occurs 
spontaneously, in response to a change in the external parameter F rather than in the structure of 
the network. 

Negative-compressibility metamaterials can be designed using this system as constituent ele- 
ments, as shown in Fig. [2}} for a square lattice configuration. We first consider this material in 
the limit of small temperatures, T, and for changes of the external force sufficiently slow that the 
lattice is in quasi-static equilibrium at any particular time (except during destabilisation transi- 
tions). Under these conditions, the material exhibits a hysteresis loop for stress-driven cycles that, 
as shown in Fig.|2j:, is identical to the one predicted for the constituents. Starting from a coupled 
configuration for all constituents (bottom left), increase in tensional stress (denoted a and defined 
as tension per lattice spacing) leads to increase in deformation, e, until a decoupling destabilisa- 
tion occurs at a = a d . A negative compressibility transition is observed at this point, where the 
material contracts by a finite amount (i.e., e decreases) upon an arbitrarily small increase of a. 
This discontinuous behaviour is indicative of a first-order phase transition in the material. The 
finite change in strain can be regarded as a "latent strain," broadly analogous to the latent heat 
absorbed during, say, the vaporization of water. Stress in our system plays the role of the temper- 
ature in this liquid- vapour analogy. The negative compressibility in our system is then the result 
of a negative latent strain, i.e., a finite drop in e upon a small increase in a. After the decoupling, 
further increase in a leads to new increase in e. More interesting, however, is the effect of reducing 
a after the transition. This eventually leads to a coupling destabilisation at a c (analogous to the 
condensation of water), which is also discontinuous, closing the hysteresis loop. For the coupling 
transition, the latent strain is positive, as a finite drop in e results from a small decrease in a, mak- 
ing this a positive compressibility transition. Each of these directional transitions is analogous to 
nonequilibrium phase transitions occurring, for example, in the Schlogl chemical reaction model 
ifTTl . and they can be interpreted as different halves of the same phase transformation, which are 
set apart by hysteresis. For an animation of the response of the material to different stress profiles, 
see Supplementary Movie. 

These destabilisations can be related to spinodals (points where metastable phases cease to 
exist), which have been considered, for example, in the stability limit hypothesis for the anomalous 



3 



behaviour of water [29 J . The possibility of destabilisations occurring on a microscopic scale and 
resulting in a phase transition has also been examined in the context of ferroelastics, such as 
for martensitic transformations and shape memory effects ll30ll . and theoretically in continuum 
mechanics PTTl and discrete systems ll32l . However, the possibility of materials exhibiting negative 
compressibility transitions like those considered here has not been previously recognised. 

A complementary but important aspect concerns the response of the material when the defor- 
mation e is taken as the tunable external parameter. This defines a hysteresis loop for strain-driven 
cycles, which is shown in Fig.[2j± Ordinarily, when the (thermodynamic conjugate variable) strain 
is taken as the controlled variable, the stress remains constant during the phase transformation as 
the latent strain is overcome through the formation of a phase mixture. In our system, such be- 
haviour does occur during the coupling transition, but during the decoupling transition, the latent 
strain is negative. In this case, once the applied strain is increased past the destabilisation at e d , the 
resulting stress a must increase discontinuously to become compatible with the new strain, giving 
rise to a mechanism for force amplification. Moreover, instead of developing a phase mixture in 
which the proportions change continuously, the only possibility for this process is that the entire 
system (or a finite fraction of it, depending on how the strain is applied) undergoes the decoupling 
transition concurrently. We note that the conclusions from Fig. [2]:,d are valid for a cycling time r 
(the time to vary the stress and strain between their minimum and maximum values) larger than 
the characteristic time scale r for the system to approach equilibrium in the limit of small T. 

Figure [3] shows how the phenomenology changes for finite temperatures. For vanishingly small 
T, the thermodynamic equilibria reduce to the configurations that are local minima of the potential, 
as considered in Fig.[2j As the temperature is increased from zero, the initially stable configuration 
around a local minimum becomes metastable. But how does hysteresis, and hence the rate at which 
the thermal fluctuations allow the system to escape a local minimum of the free energy, depend 
on the temperature? To address this question, we established a rigorous model to estimate the 
effective energy barrier for both decoupling (E%) and coupling (E d ) transitions (Supplementary 
Information). While the decay of metastability depends on the free energy, it is the height of an 
effective energy barrier — which also accounts for the frequency with which the free energy barrier 
is approached, the cycling time, and the size of the sample — that determines the transitions to the 
global minimum (Supplementary Information). For stress cycles, for example, the effective energy 
barrier for decoupling transitions in a square lattice ofNxN constituents satisfies 

2W 



EZ(a, N,r)> r^^J'"" T . <3 



log 



7r a d 



-r(u) L (a) +u T (a)) 



where V b c is a potential barrier, and uj l and uj t are the frequencies of the lattice vibration modes 
that give rise to decoupling events. To test this model, we have implemented molecular dynamics 
simulations (Methods Summary) for both stress cycles (Fig.[3^,c) and strain cycles (Fig.[3j3,d). In 
all cases, the transitions occur earlier than in the zero temperature limit. More important, the onset 
of the transitions (Fig.|3]:,d) is determined by the points at which the thermal energy k B T becomes 
comparable to the effective energy of the barrier, E^ d (Fig. 3p,b). The transitions are determined 
using a cycling time r larger than r . For nonzero T <C E% jk B , r can be interpreted as a time 
scale for the metastable state to be approached. 

The evolution of the phase transition at nonzero temperatures depends on the interplay between 
three time scales: the cycling time r, the characteristic time to approach metastable equilibrium 
t , and the characteristic time to approach thermodynamic equilibrium, r t . In general, r is larger 
than but comparable to the time scale for volume fluctuations r v , and r t is larger than but compa- 
rable to the time scale for transitions across the energy barrier. Consider, for example, the stress 
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cycles shown in Fig. [3^,c. The coupled equilibrium of the potential has lower energy for small 
stresses, while the decoupled equilibrium has lower energy for sufficiently high stresses. At some 
critical stress, the energy of the two equilibria coincide. Equilibrium statistical mechanics predicts 
that the phase transition would proceed at this critical stress. Such an equilibrium phase transition 
cannot be a negative compressibility transition, as examination of the differential change in free 
energy at fixed temperature reveals that the high-pressure phase is always also the low-volume 
phase. However, the occurrence of an equilibrium phase transition involves the assumption that 
stress is changed arbitrarily slowly. For a fixed but nonzero rate of change (characterised by r), the 
external parameter may change before the system has approached its thermodynamic equilibrium 
(characterised by r t ) or even metastable equilibrium (characterised by r ). If r is larger than but 
comparable to r , it follows from the Boltzmann distribution that the temperature-dependent time 
scale r t will be exponentially large for k B T <C E£ d 9 while r t becomes comparable to or smaller 
than r as k B T becomes comparable to E^ d . This underlies much of the properties shown in Fig. 
[3] and is expected to be representative for laboratory time scales. Larger temperatures and longer 
cycle times cause the metastable states to decay earlier in the cycle, but our molecular dynamics 
simulations show that transitions giving rise to negative compressibility and force amplification 
persist for a wide range of small temperatures and for large cycle times (Supplementary Informa- 
tion, Figs. S3 and S4). Furthermore, as indicated in Eq. (|2j), we can show that E c ^ d depends only 
logarithmically on the cycling time r and size TV of the sample, indicating that orders of magnitude 
changes in them produce only moderate changes in the effective energy barrier (Supplementary 
Information). 

This last observation is a very important. We conservatively predict that systems with typical 
atomic-sized length and energy scales can exhibit negative compressibility transitions in macro- 
scopically accessible time scales. For example, given an energy barrier typical of a covalent bond, 
AVt c ~ 1.5 eV, microscopic frequencies typical of atomic vibrations, ujl,^t ~ 10 14 s" 1 , and 
macroscopic system size, TV 2 ~ 10 23 , we find that the effective energy barrier El becomes compa- 
rable to thermal fluctuations at room temperature, k B T = 0.026 eV, only when r > 10 13 s. Thus, 
metastable states with such an energy barrier will survive thermal fluctuations for over 3 x 10 5 
years! But of course AV& C decreases in size as a approaches a d (where AVf = 0). Therefore, the 
metastable states go from very long-lived to short-lived as the stress increases, so that the decay 
of metastability produces the desired negative compressibility transition in an accessible amount 
time. Given the constituents needed to achieve negative compressibility transitions, their design 
may be more easily accomplished using mesoscopic structures, which will generally involve yet 
larger energy barriers and hence even more robust metastable states. 

Taken together, our results show that one can design metamaterials that exhibit both negative 
compressibility and force amplification transitions (Fig. [4]). The same properties can be replicated 
for different network and lattice structures and in three dimensions. Potential applications of these 
materials include the development of new actuators and protective mechanical devices, and the 
enhancement of MEMS technologies. Macroscopic implementation is achievable with established 
manufacturing techniques and could find use in applications such as the design of smart deploy- 
able space structures (33) . We also note that the mechanisms by which mechanical networks 
can expand in response to pressure and contract in response to tension admit formally equivalent 
constructions in other types of networks. We anticipate, in particular, that they can be exploited 
to create new devices with negative finite-difference resistance in electrical and microfluidic net- 
works. 
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METHODS SUMMARY 



Stable equilibrium of constituents . For a given F, the equilibrium conditions for which all forces 
are balanced are given by dV/dx = dV/dy = dV/dh = 0. If V£(x) > 0, VJ(y) > 0, and 
V{l(h) > 0, as assumed throughout, the equilibrium points (x* (F) , y* (F) , h* (F)) satisfy h* = 
x* + y*. An equilibrium point is linearly stable if the matrix of second derivatives of V, denoted 
D 2 V 9 is positive definite at that point. In addition, at stable equilibrium points dh*/dF > 0, 
confirming that negative compressibility is ruled out for infinitesimal force changes. 

Condition for destabilisation of constituents. The Sylvester's criterion applied to the Hessian 
matrix D 2 V shows that an equilibrium point is stable if and only if 

z[U ; V£{x*) + Vj(y*) + 2V»(x* + y*) K } 

(for details, see Supplementary Information). Destabilisations occur at degenerate equilibrium 
points, where the inequality is saturated. 

Potentials leading to negative compressibility. We let V z be attractive for all but small z and 
become negligible for large z 9 and we let V X9 V y , and V h be increasingly attractive for large x, y, 
and h 9 respectively. Unless otherwise noted, we considered V z to be the Lennard- Jones potential, 
Vh to be a harmonic potential, and, for conceptual simplicity, V x and V y to be piecewise harmonic 
potentials. The choice of parameters is described in the Supplementary Information, where we 
also show that similar effects can be obtained using smooth V x and V y potentials. 

Molecular dynamics simulations. The simulations were carried out on a square "AT-cell" consist- 
ing ofNxN constituents with periodic boundary conditions, and were based on a generalisation 
of the Andersen barostat ll34l (Supplementary Information). The prescribed strain is achieved by 
changing the position of the boundaries and the prescribed stress is achieved by introducing addi- 
tional degrees of freedom that account for volume fluctuations. Constant temperature is achieved 
using Nose-Hoover thermostat chains E51 . in which yet additional degrees of freedom are intro- 
duced to exchange energy with the system. 
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Negative Compressibility Transitions 
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Fig. 1. Negative compressibility contrasted with other effects, a, b, Negative compressibility transitions: 
longitudinal deformation A r L x opposing the longitudinal applied force change A a F x , which can be tension 
(a) or pressure (b), such that A r L x /A a F x < 0. Here, index a denotes "applied," index r denotes "response," 
and A a,r is positive (negative) for an increase (decrease) of the corresponding quantity, c, Negative Poisson's 
ratio: transverse expansion (contraction) A r L y in response to longitudinal stretching (shrinking) A a L x , such 
that — A r L y / A a L x < 0. d, Negative incremental stiffness: decrease in the resulting restoration force F x in 
response to an increase in deformation, A a L x , such that A r F x / A a L x < 0. This behaviour can be stable for 
one component of a composite (e.g., for inclusions in a metal matrix) but not for the material as a whole, 
e, Stretch-densified materials: materials (with positive bulk modulus) that expand in one direction when 
hydrostatically compressed also become denser when stretched along this direction, such that for small 
deformations A a L x one has (A r L y )L x + (A a L x )L y < 0. Note that the responses in (a) and (b) do not 
follow from the response in (d). The significant difference is that in (d), due to constraints necessary to 
stabilise otherwise unstable configurations, the independent variable can be taken to be the deformation 
A a L x rather than an applied force (see also Supplementary Information, Sec. 9). This distinction is crucial 
and should be contrasted with (c) and (e), where it generally makes no difference whether the applied 
quantity is the longitudinal force or the longitudinal deformation. 
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Fig. 2. Constituents, metamaterial, and resulting hysteresis loops, a, Constituents, consisting of four 
particles that interact via potentials V x , V y , V z , and Vh, can be found to be either in a coupled (red) or 
decoupled (blue) stable configuration. The realised configuration depends on the applied force and history. 
Note that z = y — x, x = h — y, and y = h — x, so that it takes just three independent variables, x, y, and 
h, to describe the system completely (see equation (1)). The system undergoes a negative compressibility 
transition as the applied force F increases and the configuration changes from coupled (left) to decoupled 
(right), b, Metamaterial consisting of a square lattice of constituents under a general application of stress, 
where the individual particles are indicated using the colour scheme of (a). There is a small but noticeable 
contraction (of w 5%) at the centre of the sample, where the tension is higher (the force gradient is exag- 
gerated to facilitate visualisation). The white background marks the boundaries of the material prior to the 
contraction, c, d, Hysteresis loops of the constituents (dotted green) and metamaterial (continuous purple) 
for tunable uniform applied tensional stress a (c), and for tunable uniform applied strain (deformation) e (d). 
For vanishingly small temperature, as considered in this example, the loops of the material and constituents 
differ only at the coupling destabilisation of strain cycles, which is discontinuous for each constituent and 
continuous for the material in the thermodynamic limit. In (c) and (d), the small displacements between 
the curves were introduced to facilitate visualisation. Because the equilibria of the constituents involve 
two rather than one degree of freedom, the crossing points of the loops correspond to different equilibria 
configurations, one coupled and the other decoupled (Supplementary Information, Fig. SI). 
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Fig. 3. Effective energy barrier and the temperature dependence of the transitions, a, b, Effective 
energy barrier between the coupled and decoupled configurations of the constituents, E^ d , as a function 
of the applied stress (a) and applied strain (b). The height of the barrier is measured relative to the stable 
coupled state (dotted red; E%) and to the stable decoupled state (continuous blue; Ejf). c, d, Fraction of 
the material in the coupled phase, x> f° r several temperatures as a function of the applied stress (c) and 
applied strain (d). The curves were generated by averaging over 10 realisations of our molecular dynamics 
simulations on a square lattice of 20 x 20 constituents with constant rate of change of the tunable parameter 
(a or e) for a cycling time r = 10 L5 tv, where ry is the time scale of volume fluctuations. The small vertical 
displacements in the curves in (c) and (d) were were introduced to facilitate visualisation. The temperature 
To is a reference, for which fcgTo is the potential barrier separating the coupled and decoupled phases when 
they have the same potential energy. The horizontal dotted lines indicate the points where the energy barrier 
equals the thermal energy fcgT. It is observed that these points are excellent predictors of the onset of all 
coupling and decoupling transitions in (c) and (d) (vertical dotted lines). For details on the effective energy 
barrier, reference temperature, and volume fluctuation time scale, see Supplementary Information. 
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Fig. 4. Physical examples of stress-induced negative compressibility and strain-induced force ampli- 
fication transitions, a, A rope made of our material holds a basket containing a weight, which stretches 
the rope to near the decoupling transition point. When a small amount of extra weight dW is added to the 
basket, an instability develops, and the rope spontaneously passes from the coupled to the decoupled phase. 
As a result, it raises the basket, despite the increase in the weight, b, A rope made of our material is attached 
to the ground and to a raised ratchet. The rope is stretched to a length near the decoupling transition point. 
When the ratchet is rotated a small amount d9, an instability causes the rope to pass from the coupled to the 
decoupled phase. As a result, the tension in the rope is amplified by a large amount. 
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